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Abstract 

The relevance of the quasi-two-dimensional spin-1/2 frustrated quantum antifer¬ 
romagnet due to its possibility of modelling the high-temperature superconduct¬ 
ing parent compounds has resulted in numerous theoretical and experimental 
studies. This paper presents a detailed research of the influence of the varying 
exchange interactions on the model magnetic properties within the framework of 
self-consistent spin-wave theory based on Dyson-Maleev representation. Beside 
the nearest neighbour interaction within the plane, the planar frustration up to 
the third nearest neighbours, cyclic interaction and the interlayer coupling are 
taken into account. The detailed description of the elementary spin excitations, 
staggered magnetization, spin-wave velocity renormalization factor and ground- 
state energy is given. The results are compared to the predictions of the linear 
spin-wave theory and when possible also to the second-order perturbative spin- 
wave expansion results. Finally, having at our disposal improved experimental 
results for the in-plane spin-wave dispersion in high-Tc copper oxide La 2 Cu 04 , 
the self-consistent spin-wave theory is applied to that compound in order to cor¬ 
rect earlier obtained set of exchange parameters and high temperature spin-wave 
dispersion. 
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interaction, Quasi-two-dimensional antiferromagnetic copper oxides 


1. Introduction 


The importance of the quasi-two-dimensional (2D) antiferromagnetically 
(AFM) ordered systems, corroborated by the fact that the high-temperature 
superconducting cuprate parent compounds, La 2 Cu 04 for instance, belong to 
that class of systems, highly justifies their decades-long presence in the scientihc 
literature. The interest is especially devoted to their magnetic properties in or¬ 
der to elucidate the underlying microscopic mechanism of the high-temperature 
superconductivity. Though the 2D spin-1/2 Heisenberg model with only the 
nearest neighbour (NN) interaction coupling J was frequently used as the re¬ 
search initial point [l|, numerous studies, both theoretical and experimental, 
have in the meanwhile indicated the significance of taking into account other 
exchange interactions (see [ 3-01 and references therein), leading to the study 
of the frustrated spin systems |l3|. Most emphasized is the relevance of the 
planar frustration, induced by the next-(NNN) or next-next-nearest (NNNN) 
neighbour interaction J 2 / 3 , as well as the cyclic exchange Jc comprising four- 
spin plaquette interaction, often recognized as the second strongest exchange 
interaction in the system. 

The self-consistent spin-wave theory (SCSWT) with its numerous versions is 
widely used for the analysis of the Heisenberg layered magnets (see [lll4l4l| and 


references therein). Our research was primarily motivated by the study of the 
three-dimensional (3D) antiferromagnet with the planar frustration and cyclic 
exchange interaction within the SCSWT performed in Ref. 0 . The theory 
was applied strictly to La 2 Cu 04 , without giving a detailed study of the depen¬ 
dence of the model magnetic properties on the varying values of the exchange 
interactions contained in Hamiltonian. This, however, may be interesting to 
investigate, especially since it is experimentally possible to vary the frustration 
ratio J 2 /J by applying high pressures on the system 1 ^. The present paper 


is intended to be a comprehensive study on the influence of planar frustration 
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ratio, cyclic exchange and interplanar exchange interaction on the spin-wave 
spectrum, ground-state sublattice magnetization, spin-wave velocity renormal¬ 
ization factor and ground-state energy of the three-dimensional (3D) quantum 
antiferromagnet with S' = 1/2 within the self-consistent spin-wave theory. We 
shall determine the exchange parameter region where the ground state of the 
Neel type exists. The existence of the other long-range ordering patterns (colum¬ 
nar phase, for instance), as well as the intermediate quantum disordered phases, 
will not be the issue of this paper. The SCSW approach will be compared to the 
linear spin-wave theory (LSWT) result^ as well as the results of second-order 
perturbative spin-wave expansion from [6l-l9| . The influence of the third nearest 
neighbour interaction on sublattice magnetization will be analyzed, which to 
our knowledge has not been earlier discussed within the framework of SCSWT 
applied on the above introduced model. Special attention will be devoted to the 
high-Tc superconducting parent compound La 2 Cu 04 , thoroughly studied in the 
literature. Using the improved experimental in-plane spin-wave dispersion 
we shall find the corrected exchange parameter set, whereby we shall calculate 
spin-wave dispersion at high temperatures and compare it both to the experi¬ 
ment from the Ref. Q], as well as to the numerically obtained spectrum from 

El- 

The paper is organized as follows: in Section 2 we present the model Hamil¬ 
tonian and introduce the dominant exchange interactions. In Section 3 the 
Hamiltonian is mapped into the bosonic form based on the Dyson-Maleev (DM) 
transformation. Mean-field decoupling of the fourth-order terms is performed 
and main expressions for the quantities to be analyzed are derived by making 
use of the self-consistent spin-wave theory. The numerically calculated results 
are plotted and analyzed in detail in Sec. 4. The conclusions are shortly stated 
in Sec. 5. 


3 



2. Model Hamiltonian 


The dominant interactions for the S' = 1/2 antiferromagnet on the tetragonal 
lattice are presented in Figure 1 and comprised in the following Heisenberg 
Hamiltonian: 




1 


+5i 


JiSf, 


■ ^fia.+Si + 


na ,S-l_ 


(^ = a.b) 


(i=2,3) 


+Jc .( 

( 


Sn. ■ S, 


fla -\-di 


Q ^ . Q 

Ticx-\-d2 n.Q.+d 3 


)-( 


s«„ • s. 


rZc+da 


)( 


^'fT-a~\~di ^f^a~\~d2 j 


(c. = a,b) 

- { 


na+d 2 


)( 


5 




)]+j^ E E 


•s 


na-\- 6 ‘^ 




( 1 ) 


The first two terms include the interactions between the first, second and third 
nearest neighbours within the plane. The position of the spin in the sublattice 
a (a = a, b) is denoted by ria, while the vectors connecting the given spin and 
its corresponding neighbours are denoted by 5i {i = 1,2,3). The third term 
presents the cyclic exchange interaction, described by the exchange integral 
Jc- The vectors connecting the spins in a plaquette read: di = a (e*r + Cy), 
d 2 = a Cx and ds = a (e^ — Sy)- The interplanar interaction is described by the 
last two terms, where the vectors designate the vectors connecting the 

antiferro-/ferromagnetically coupled nearest neighbour spins in adjacent planes. 
All interactions are assumed to be antiferromagnetic, i.e. J, J 2 , J 3 , Jc, J± > 0. 
Contrary to Refs, [s-^, we observe that Hamiltonian (HD includes the next-nex' 
nearest neighbour planar interaction, enabling us to incorporate in the present 
study the influence of this interaction in the Neel phase. 
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Figure 1: (Color online) Magnetic unit cell (bold solid line) of 3D antiferromagnet 
with exchange interactions labeled. Two different orientations of spins are denoted 

by open and solid circles. 

It is common to define the fundamental energy scale by the NN exchange 
interaction J. Therefore, we introduce the dimensionless ratios and study their 
influence on the spin-wave spectrum, ground-state sublattice magnetization etc. 
These quantities are the frustration ratio A 2/3 = JiizlJ-, the ratio Ac = Jc/J 
parametrizing the relative strength of the cyclic exchange and the ratio A_l = 
Jj_/J parametrizing the relative strength of the interplanar interaction. 


3. Spin-wave spectrum, magnetization and related quantities 


In order to determine the spin-wave spectrum, we apply a variant of the 
non-linear spin-wave theory [l 6 |. The spin operators in Hamiltonian ([IJ are 


first written in the local coordinate system 


17| and then expressed as boson 


operators using the well-known Dyson-Maleev transformation in the form 
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In terms of these operators, the Hamiltonian ([1]) is mapped into an equivalent 
Hamiltonian of interacting bosons, which may be written as 


H = Ho+H 2 +Hi, 


(3) 
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where Hq presents the classical ground-state energy, H2 denotes the bilinear 
part of the Hamiltonian which corresponds to the non-interacting spin waves 
(LSW), while describes the interaction among spin waves induced via trans¬ 
formations (m and consists of the terms quartic in Bose operators. The terms of 
the higher order are immediately neglected. We then simplify the Hamiltonian 
© in the spirit of the Wick’s theorem in the manner similar to the one 
implemented in Ref. . The manner in which the quartic terms are evaluated 
is illustrated on the following example: 
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where the single brackets (...) indicate averages over the canonical ensemble at 
temperature T. As a consequence of DM transformation, the obtained Hamil¬ 
tonian is non-Hermitian. Therefore, we symmetrize it by adding the Hermitian 
conjugate terms to the non-Hermitian ones [^. After decoupling procedure is 
completed, we perform the Fourier transform and obtain the following bosonic 
Hamiltonian 


H = -JS^z— + J 2 S^—z + J 3 S^—z + J^S^N + 

+ 

k k 

Here, z denotes the coordination number and equals four in all cases, N presents 
the number of lattice sites, while the quantities e{k) and a{k) are given by the 
expressions 


z{k) = JSz Ti - A2r2(l - 72(fc)) - A2r3(l - 73(fc))- 


- + r,( 2 ) 72 (fc)) + Ai(rf - rr(i - 77 (fc))) 


Y(k) = JSz 


ri7i(fc) - yr,( 3 ) 7 i(fc) + A^rf 7 f (fc) 


( 6 ) 


(7) 


The quantities 7(fc) denote the structure factors with lattice constants set to 
unity: 

k k ^ 

7i(fc) = cosycos-^ 72(fc) =-(cos fca;-I-cos fcy) 'y3{k) = coskxcosky , 
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while the quantities F present the renormalization factors given by 
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(9) 


In Bloch’s approximation all these renormalization factors equal unity, while 
here due to the quantum and thermal fluctuations they are defined by the afore¬ 
mentioned expressions. The temperature functions which enter Eqs. ([9|) are the 
following correlation functions 
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where 9 = kgT. The diagonalization of ([5]) by making use of the standard 
Bogoliubov’s transformation yields the spin-wave spectrum 

E{k) = \Je‘^{k)-a^ik), (11) 


while the ground-state energy per lattice site, with the quadratic and quartic 
corrections included, reads 


Eo/N = -^JS^z + ^J 2 S^z + ^sS^z + ^ ■ ( 12 ) 

k 


The spectrum given by m possesses the Goldstone mode, as expected due to 
the spin isotropy of the model. 

The sublattice magnetization in the absence of the external magnetic field 
is given by 


Na^ 

k 


iiW COthSi _ i 

2 E{k) 20 2 


(13) 


Finally, we shall examine the spin-wave velocity renormalization factor Zc 
in 2D case. Namely, in the long-wave limit, the spin-wave dispersion vanishes 
linearly according to 


E(k) ^ ck, k = ^Jk"^ + ky —>■ 0 , 


(14) 


where c presents the spin-wave velocity given by 


c = 2V2ZcJS. 


(15) 


Therefrom, the renormalization factor Zc reads 

Zc = ^J [hi - _ 2A2r2 - dAaTg + 2A,52(r^(2) - r.^))]. 

(16) 

All the numerical calculations will be performed based on upper expressions. 
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4. Analysis of results 


4.I. Spin-wave energies 

In order to describe the ground-state behavior of the system we perform the 
calculation of the spin-wave spectrum for the 2D model within LSW and SCSW 
theory, based on numerical evaluation of Eqs. ®, © and (EH), where we put 
T = 0. In the LSW approach, T quantities in last two expressions equal unity, 
while within the SCSW approach the cumbersome system of Eqs. ®-Eni) has to 
be solved. Due to its self-consistency, the iterative procedure has to be applied. 
In Eigure 2 we show the comparison between the spin-wave spectra along the 
high symmetry directions in the 2D AFM Brillouin zone (for the tetragonal 
phase) obtained within the LSW (blue line) and SCSW (red solid line) theory, 
for the relative frustration A 2 = A 3 = 0.1 and cyclic exchange Ac = 0 . 



Figure 2: (Color online) Spin-wave energy E{k)/J along the high symmetry 
directions in the 2D AFM Brillouin zone obtained from LSWT, SCSWT and first 
correction to LSWT (red dashed line), for the couplings A 2 = 0.1 and Ac = 0. 


We also show the results obtained as the first correction to the LSWT result, 
i.e. by inserting the LSW results in the SCSW expressions, without performing 
the self-consistent procedure (presented in Fig. 1 with the red dashed line). If 


B3 


not stated otherwise, hereafter we assume that A 3 equals A 2 lau and do not 
emphasize the value of the former for brevity. While qualitatively LSWT and 
SCSWT results resemble each other, the SCSWT gives significant contribution 
to the entire dispersion curve. Further, it can be seen that the SCSWT results 
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differ only slightly from the first correction to the LSWT and only in the vicinity 
of the Brillouin zone boundary. Other choices of exchange parameters confirm 
this observation, whereby this difference grows with the increase of Ac, though 
very slowly, being less than 1% for A 2 = 0.1 and Ac = 0.6, which may bring us to 
the conclusion that the tedious procedure of solving the aforementioned system 
of equations iteratively is not justified. Though, the study of the ground-state 
sublattice magnetization and other quantities of interest will show that it is 
necessary to employ the SCSW theory in order to avoid the divergences which 
arise from taking the first correction to the LSWT only. 

In Figure 3a) we present the LSW theory results for the spin-wave spectra 
for different values of parameters A 2 and Ac. 




a) b) 

Figure 3: (Color online) a)Spin-wave energies E{k)/J obtained within LSWT, for 
various values of parameters A 2 and Ac. b)Spin-wave energies E{k)/J within 
SCSWT, for various values of parameters A 2 and Ac. 

It is obvious that with the increase of frustration, the magnon energies get 
lowered, with a clearly observable dip from (tt, —tt) to (tt, 0 ), meaning that the 
excitation energy at (tt, 0) decreases more strongly than at (tt, —tt). The dip 
grows with the frustration and is absent only for A 2 = 0 and Ac = 0, as also 
stated in Ref. Q. The increase of cyclic exchange for fixed A 2 supports the 
softening of (tt, 0 ) mode, while the (tt, —tt) mode remains unaffected, which 
presents the shortcoming of the LSW approach. This drawback is removed 
within the SCSWT, as seen in Figure 3b), where the evolution of the spin-wave 
spectra with varying parameters A 2 and Ac is plotted. Within SCSWT the 
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increase in Ac softens also the (tt, —tt) mode, in accordance with [7|. We observe 
that even in the absence of the cyclic interaction there exists a dip from (tt, —tt) 
to (tt, 0) (except for A 2 = 0), which grows with frustration. Hence, in this case, 
SCSW theory gives result opposite to the one obtained within the second-order 
perturbative spin expansion , which predicts the dip of the same sign only for 
the substantial cyclic exchange, while for Ac = 0 the inverse dip from (tt, 0 ) to 
(tt, —tt) is obtained. 

4-2. Ground-state sublattice magnetization 

We proceed with the calculation of the ground-state sublattice magnetiza¬ 
tion, again for the 2D model at first, by making use of the expression m, 
whereby we again have to apply the iterative procedure. The LSWT and SC- 
SWT results for the zero-temperature magnetization dependence on the frus¬ 
tration ratio are compared in Figure 4a). 




a) b) 

Figure 4: (Color online) a)The ground state sublattice magnetization dependence on 
A 2 for different values of Ac. The solid lines present the LSWT and SCSWT results. 

The dashed lines present the results of the first correction to the LSWT. b)The 
ground state sublattice magnetization vs. Ac for different values of A 2 in LSW and 

SCSW approach. 

Since the presence of frustration disorders the system, the sublattice magnetiza¬ 
tion {S^)o is gradually reduced by the increase of A 2 , i.e. the zero-point quantum 
fluctuations grow. The magnetization within SCSWT is enhanced and decreases 
more slowly compared to the LSWT predictions, in accordance with Q. The 
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apparent divergence which arises within the non-linear SW theory including the 
first correction to I/S' order only (shown in Figure 4a) by the dashed lines), 
is successfully removed by allowing for the SCSW approach. This conclusion 
agrees with those quoted in The sublattice magnetization dependence 

on the cyclic exchange parameter is presented in Figure 4b). It can be seen 
that the cyclic exchange stabilizes the system, whereby the magnetization is 
enhanced within SCSWT. Besides, the difference between the LSW and SCSW 
theory results grows with the increase of the parameter 

A2. 


It is important to emphasize that the domain of the parameters A 2 and 
Ac where the ground state of the Neel type exists for the model defined by 
the Hamiltonian o within SCSWT, is limited. The corresponding parameter 
region is shown in the phase diagram in Figure 5. 



Figure 5; (Color online) The Neel phase parameter region in the J-J 2 -J 3 -JC model. 


Further, we study the influence of the third nearest neighbour interaction on 
the sublattice magnetization, in order to estimate its contribution to the ground- 
state sublattice magnetization. Though the NNNN interaction was included in 
the Hamiltonian in Ref. its influence on the model magnetic properties 
was not separately discussed. The comparison between the LSWT and SCSWT 
results for the cases with and without the NNNN interaction (for simplicity in 
the absence of cyclic exchange) are presented in Figure 6a). 
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Figure 6: (Color online) a)The comparison of ground-state sublattice magnetization 
dependence upon A 2 for A 3 = 0 (dotted lines) and A 3 = A 2 (solid lines). For 
simplicity, Ac = 0. b)The ground-state sublattice magnetization dependence on the 
parameter Ax, for two choices of parameters A 2 and Ac. 


Evidently, with the presence of the non-vanishing parameter A 3 the frustration in 
the antiferromagnet grows, additionally destabilizing the system. For instance, 
within the LSWT, the critical value of A 2 at which the AFM order is destroyed, 
equals approximately 0.38 in the absence of NNNN interaction, which presents 


the result also stated in [ 6 |, while it falls to a much lower value of 0.13 for A 3 0. 

The SCSWT results also corroborate this observation. 

Next, we allow for the interplanar interaction between the adjacent planes 
which complicates the analysis significantly, since it introduces the additional 
frustration into the system (see Figure 1). In order to examine the effects 
of the three-dimensionality on the behavior of the system, we determine the 
influence of the parameter Ax on the ground-state sublattice magnetization. The 
calculated results are plotted in Figure 6 b) both for the simple model with A 2 = 
A3 = Ac = 0 and the model with non-vanishing planar frustration and cyclic 
exchange. It may be seen that LSW and SCSW theory show opposite tendency 
to each other. Though LSWT predicts the decrease in the magnetization with 
the grow of the parameter Ax due to the frustration which destabilizes the 
system, SCSW theory gives the growth in magnetization, as a consequence 
of the subtle interplay of the competing interactions. Besides, we notice that 
the same change in the parameters A 2 and Ac yields within LSWT only a slight 
decrease of magnetization, while the increase within SCSWT is more significant. 
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4-3. Other quantities 

Further, we calculate the spin-wave velocity renormalization factor by mak¬ 
ing use of expression (Hg. We first increase the frustration ratio A 2 and calculate 
the renormalization factor for different cyclic exchange ratios Ac in both LSWT 
and SCSWT approach (Figure 7a)). 



Figure 7: (Color online) a)Renormalization factor Zc for the spin-wave velocities vs. 
frustration A 2 for different values of cyclic exchange Ac, within LSWT and SCSWT 
approach. The dashed lines present the results of the first correction to LSWT. 
b)Renormalization factor Zc vs. Ac, for different A 2 . 

The results demonstrate that Zc decreases with A 2 , the decrease being less 
rapid within the SCSW theory. The divergence similar to the one appearing in 
the analysis of (<5^)o is again eliminated within SCSWT. The renormalization 
factor is also calculated in dependence of Ac for different values of A 2 (Figure 
7b)), showing the more rapid decrease of Zc within SCSWT. 

Finally, by making use of the expression (USD, we analyze the ground-state 
energy per lattice site versus frustration and cyclic exchange parameters. The 
dependence on the parameter A 2 is plotted in Figure 8a). 
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Figure 8:(Color online) a)Ground-state energy per lattice site vs. A 2 , for different Ac, 
in LSW and SCSW approach. b)Ground-state energy per lattice site vs. Ac, for 

different A 2 . 


It is interesting to notice that the absolute value of the ground-state energy 
within SCSWT at first decreases with frustration, however, after a certain value 
of the frustration ratio it starts to increase. Figure 8b) presents the dependence 
of the ground-state energy on the parameter Ac. The results show that \Eo/N\ 
decreases with the growth of the cyclic exchange parameter. 


2 ^.^. Application to La 2 Cu 04 

The model described by the Hamiltonian (ED may be used to study the 
magnetic properties of the high-temperature superconductor parent compound 
La 2 Cu 04 . This compound has been the subject of our previous studies [l^, 
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22|, where the calculation was performed within the method of spin Green 


functions on the model without the cyclic exchange interaction. The model with 
the non-vanishing parameter Ac within the SCSW theory has been examined in 
detail in Ref. The set of bare superexchange couplings quoted therein 

was determined by the accurate fit of the in-plane spin-wave dispersion from 


Ref. 0, obtained by the inelastic neutron scattering at 10 K. Since we had 
improved experimental data from Ref. Q at our disposal, we recalculated the 
exchange parameters. The frustration parameter is taken to be as small as 


A 2 = 0.01, according to Q, instead of the value 0.025 from [^. The set of the 


a. 


exchange parameters is thus slightly changed and within the SCSW theory reads: 
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J = 142.62 meV, Ac = 0.22 and A 2 = A 3 = 0.01. If we use this set of parameters 
to obtain the spin-wave dispersion at 295 K, we obtain the spectrum which 
gives up to 25% of the observed changes between the low- and high-temperature 
spectrum. This presents a mild improvement to the spectrum calculated from 
the parameters quoted in since the later accounts only for a few percent of 
the observed changes. 


5. Conclusions 

The three-dimensional tetragonal S' = 1/2 antiferromagnet presents an in¬ 
teresting system, studied in literature due to its application to layered copper 
oxides, which present the high-temperature superconducting parent compounds. 
In the present paper this system is described by the Heisenberg Hamiltonian 
comprising the planar interactions up to NNNN interaction, the cyclic exchange 
interaction, as well as the interplane interaction. The magnetic properties of the 
ground state of the model are examined within the framework of SCSW theory 
and compared to the LSWT predictions. 

We show that the system is highly frustrated and therefore subjected to the 
complex competition of different exchange interactions, leading to the specific 
magnetic properties behavior. The spin-wave spectra depend on frustration and 
cyclic exchange ratio, whereby the increase of these ratios yields the softening 
of the spectra. The sublattice magnetization dependence on the aforementioned 
parameters demonstrates that the zero-temperature quantum fluctuations play 
a pronounced role. The planar frustration and cyclic exchange have the opposite 
impact on the ground-state sublattice magnetization, whereby the former desta¬ 
bilizes, while the later stabilizes the system. The effect of the next-next-nearest 
neighbours interaction is shown to be significant, according to the fact that the 
critical value of the frustration ratio where AFM order is destroyed is, within 
LSWT for instance, three times smaller in the presence of the NNNN interac¬ 
tion. The influence of the interlayer coupling on the sublattice magnetization 
is also examined and turned out to be dictated by the exchange interaction 
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competition. The spin-wave velocity renormalization factor and ground-state 
energy are strongly affected by the frustration and cyclic exchange interaction 
as well. 

The SCSW theory is also applied to layered cuprate oxide La 2 Cu 04 . The set 
of exchange parameters in this compound is determined to fit the experimental 
spin-wave dispersion obtained at T = lOK [^. The high-temperature spectrum 
is then calculated, showing better agreement with the experiment than the 
earlier ones, but still indicating to the limitations of the spin-wave theories at 
high temperatures. 

The aim of this study was to give an exhaustive research of the results which 
SCSW theory gives for this model, as well as to point out its shortcomings, 
where some more sophisticated theories, as for example perturbative spin-wave 
expansion up to 1/S'^, have to be applied. 
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